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Customization for Cryogenic Tank Analysis «***» 


• A general purpose SINDA/FLUINT (S/F) stratified tank model was 
created to simulate self-pressurization and axial jet TVS 

• Stratified layers in the vapor and liquid are modeled using S/F 

lumps. 

• The stratified tank model was constructed to permit incorporating 
the following additional features: 

- Multiple or singular lumps in the liquid and vapor regions of the 
tank 

- Real gases (also mixtures) and compressible liquids 

- Venting, pressurizing, and draining 

- Condensation and evaporation/boiling 

- Wall heat transfer 

- Elliptical, cylindrical, and spherical tank geometries 
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Customization for Cryogenic Tank Analysis «***» 


• Extensive user logic is used to allow detailed tailoring - Don't 
have to rebuilt everything from scratch!! 

• Most code input for a specific case is done through the 

Registers Data Block: 

- Lump volumes are determined through user input: 

• Geometric tank dimensions (height, width, etc) 

• Liquid level could be input as either a volume 
percentage of fill level or actual liquid level height 





SINDA/FLUINT Stratified Tank Setup 

• • * l 


The number of lumps in the vapor 
1 region need not be the same as the 
number of lumps in the liquid region. 


The number of wall nodes in the 
2 vapor region need not be the same 
as the number of wall nodes in the 
liquid region. 



Number of Wall Vapor Nodes = 3 


Number of Wall 
Liquid Nodes = 3 



Number of Vapor 
Lumps = 3 (Equal 
Volume 


S/F I FACES 
between ALL 
lumps 


UQUIO LEVEL 


Liquid level input as volume 
4 percentage of fill level or actual 
liquid level height. 


Number of Liquid 
1 Lumps = 4 (Equal 
Volume) 


Fig 1: Stratified Tank Setup 


LUMP VOLUMES: 

• Initial values (guesses) required 

• Extensive logic determines 
equal volume distribution 
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• Initialization of wall node volumes: 

- Determined by the number of nodes input for the vapor 
wall and liquid wall regions 

- Wall node volumes are not uniform values in the vapor 
region or liquid region 

- Wall node volumes correspond to "equal" fluid volumes in 
the vapor or liquid region 


SINDA/FLUINT Stratified Tank Setu 
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SINDA/FLUINT Stratified Tank Setup 


• Volume Flow Rate Connectors (VFRSETs) are placed between the 
layered lumps in each region to equalize volumes: 

- Used during venting, draining 

- Used during mass transfer of evaporation, condensation, boiling 

- User logic determines the volume flow rate through these 
connectors so that the volumes of the lumps in each respective 
region remains constant within a small percentage 

- Avoids lumps from becoming too small or too large relative to 
one another, and consequently avoids the problem of having 
very small lumps to somehow "disappear" 

- Converted "balloon-like-behaving" lumps with IFACES to 
something more like a "fixed volume" approach (still at constant 
pressure) 
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Wall Heat Flux (Can Be Nonuniform) 



Fig 2: Stratified Tank Heat Transfer 
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IMPOSED 
HEAT FLUX 



TANK WALL THERMAL NODE 



FREE CONVECTION HEAT FLUX TO 
FLUID LUMP 


AXIAL CONDUCTOR 


EVAPORATION AND 

LIQUID LEVEL CONDENSATION MASS 
TKANSfER INCLUDED 


S/F QTIES (2 Sets): 

• "Mixing Heat Transfer" Nu ~ Ra 0 33 

- Based on the AT Between Lumps 

• Boundary Layer mdot*Cp*AT 

- One Directional -> "UP" 


Fig 3: Stratified Tank Heat and Mass Transfer 
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• Thermal boundary layer that forms along the tank wall, due 
to wall heat leak, is modeled empirically using correlations for 
free convection. 

* Although the two dimensional flow dynamics of the 
boundary layer could not directly be incorporated into the 
one dimensional stratified tank model, it is incorporated in a 
one dimensional sense. Within each fluid lump the following 
boundary layer characteristics are determined: 

- Characteristic velocity 

- Boundary layer thickness 

- Buoyancy driven volume flow rate 
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• SINDA/FLUINT contains pre-built utility functions to model 
heat and mass transfer between a liquid and vapor interface 

(TWIN TANKS) 

• For more modelling flexibility, user logic could incorporate the 
necessary physics to model a wide variety of scenarios: 



SINDA/FLUINT Vapor Liquid Interface Modelin 



LIQUID LEVEL -> VAPOR/LIQUID INTERFACE: 

• Modelled as a PLENUM (Boundary State) 

• TSAT = Saturation Temperature of "Top" 
Liquid Lump (Updated Every Iteration) 




A S/F FTIE is placed between the "bottom 
vapor lump (Tvap) and this boundary 
PLENUM 

A S/F FTIE is placed between the "top" 
liquid lump and the boundary PLENUM 


Fig 4: Vapor Liquid Interface Model 
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SINDA/FLUINT Vapor Liquid Interface Modelin 



The heat rates for these FTIEs are defined as follows: 


QFTIEvap 

QFTIEIiq 

where, 

h 

Nu 

const 


- h * A * (t 

11 Interface \ 1 

= h * A * (T 

11 m INTERFACE V 1 


VAP “ "^INTERFACE ) 
INTERFACE - "^Uq) 


- Nu / D| NTERFACE * k 
= const* Ra (1/3) 

~ 0.04 * Function(Height/Diameter) 


The net evaporation rate at the interface is calculated to be: 
mdot EVAP = (QFTIEvap- QFTIEIiq) / (heat of vaporization) 
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SINDA/FLUINT Vapor Liquid Interface Modeling 


• Evaporation is modeled via the process: 

— Liquid leaves the "top" liquid lump and enters a DUMMY PLENUM 

- Vapor enters the "bottom" vapor lump from a DUMMY vapor PLENUM 

• Condensation is modeled via the process: 

— Vapor leaves the "bottom" vapor lump and enters a DUMMY PLENUM 

- Liquid enters the "top" liquid lump from a DUMMY liquid PLENUM 

• The DUMMY PLENUMS are set at the saturation temperature of the "top" 
liquid lump (updated every iteration) 

NOTE! 

• The mass flow rate, whether condensing or evaporating, should be at 
saturated conditions. However the "bottom" vapor and "top" liquid 
lumps may not be saturated. Thus when mass is removed from either of 
these lumps the amount of energy leaving these lumps needs to be 
adjusted to account for this discrepancy. 
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Validation Cases 



* K-Site LH2 Self-pressurization (lg) 
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K-Site LH2 Self-Pressurization (lg) 




1. Flightweight insulated aluminum ellipsoidal tank 

• Internal volume: 175 ft 3 

• Tests conducted in vacuum chamber. 

• Tank is supported by 12 fiberglass composite struts. 

• Test article is enclosed by a cryoshroud whose 
temperatures are maintained with electrical heaters. 

• Tank insulated with 2 blankets of MU. 

2. Test fluid is liquid hydrogen 

3. Steady boil-off test and measurement performed at 95% fill 
and 117 kPa. 

4. Tank fill level was reduced to desired fill level. 

5. Several hours of additional venting at 103 kPa to achieve 
stationary state. 

6. Self-pressurization tests were initiated from stationary 
stratified state. 

NASA TM-103804, 1991 and NASA TM-105411, 1992 
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PRESURE (PSI) 


K-site LH2 lg Self-pressurization: zxperirnental and 
SIN DA/FLU I NT RESULTS 


SINDA/FLUINT: 75 LIQUID LUMPS, 40 VAPOR LUMPS, 50 LIQUID WALL NODES, 40 VAPOR WALL 
nodes Fluent = lumped ullage model 


KSITE 29% FILL LEVEL 
HEAT LEAK 30.0 W 



TEST DATA 

SINDA/FLUINT 

FLUENT 


ID 


PRESURE (PSI) 


K-site LH2 lg Self-pressurization: zxperirnental and 
SIN DA/FLU I NT RESULTS 


SINDA/FLUINT: 75 LIQUID LUMPS, 40 VAPOR LUMPS, 50 LIQUID WALL NODES, 40 VAPOR WALL 
nodes Fluent = lumped ullage model 


KSITE 49% FILL LEVEL 
HEAT LEAK 30.0 W 



TEST DATA 

SINDA/FLUINT 

FLUENT 


PRESURE (PSI) 


K-site LH2 lg Self-pressurization: zxperitnental and . 
SINDA/FLUINT RESULTS 


SINDA/FLUINT: 75 LIQUID LUMPS, 40 VAPOR LUMPS, 50 LIQUID WALL NODES, 40 VAPOR WALL 
nodes Fluent = lumped ullage model 


KSITE 83% FILL LEVEL 
HEAT LEAK 30,0 W 



TEST DATA 

SINDA/FLUINT 

FLUENT 


I ( 


PRESURE (PSI) 


K-site LH2 lg Self-pressurization: zxperirnental and 
SIN DA/FLU I NT RESULTS 


SINDA/FLUINT: 75 LIQUID LUMPS, 40 VAPOR LUMPS, 50 LIQUID WALL NODES, 40 VAPOR WALL 
NODES 


KSITE 29% FILL LEVEL 
HEAT LEAK 49.35 W 



TIME (HR) 


TEST DATA 

SINDA/FLUINT 


O 


PRESURE (PSI) 


K-site LH2 lg Self-pressurization: zxperitnental and . 
SINDA/FLUINT RESULTS 


SINDA/FLUINT: 75 LIQUID LUMPS, 40 VAPOR LUMPS, 50 LIQUID WALL NODES, 40 VAPOR WALL 
NODES 


KSITE 49% FILL LEVEL 



1 23456789 10 11 12 

TIME (HR) 


TEST DATA 

SINOA/FLLUINT 
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PRESURF (P5I) 


K-site LH2 lg Self-pressurization: zxperirnental and 
SIN DA/FLU I NT RESULTS 


SINDA/FLUINT: 75 LIQUID LUMPS, 40 VAPOR LUMPS, 50 LIQUID WALL NODES, 40 VAPOR WALL 
NODES 


KSITE 83% FILL LEVEL 
HEAT LEAK 49.35 W 



TEST DATA 

SINDA/FLUINT 
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PRESURE (P SI) 


K-site LH2 lg Self-pressurization: zxperirnental and 
SIN DA/FLU I NT RESULTS 


SINDA/FLUINT used 50 LIQUID LUMPS, 40 VAPOR LUMPS 

49% Liquid Fill Level. Total Heat into Tank = 30 W Fluent = lum P ed ulla S e model 


KSITE 49% FILL LEVEL 
HEAT LEAK 30.0 W 



S/F 50L40V 

— TEST DATA 
FLUENT 
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PRESURE (PSI) 


K-site LH2 lg Self-pressurization: zxperitnental and . 
SINDA/FLUINT RESULTS 


SINDA/FLUINT used 50 LIQUID LUMPS, 40 VAPOR LUMPS 
49% Liquid Fill Level. Total Heat into Tank = 49.35 W 


KSITE 49% FILL LEVEL 
HEAT LEAK 49.35 W 



S/F 50L40V 

TEST DATA 
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Validation Cases BACK UP SLIDES!! 




• K-Site LH2 Axial Jet (lg) 
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K-site LH2 Axial Jet Experiments (lg) 


mi'ASA 


Same Tank as K-site LH2 lg Self-pressurization experiments. 
Pump and jet nozzle (mixer unit) was hardware designed for 
Shuttle Centaur LH2 tank and installed in K-site LH2 tank. 


Jet nozzle and location not changed during axial jet runs, but 
the jet flow rate was varied. 


Only considering test runs where self-pressurization was used 
to pressurize tank before turning on jet (Test Series A and B). 
Tank typically pressurized to 186 kPa before initiating jet. 


Experimental data is 
available for: 
tank heat load, 
ullage pressure, 
fluid temperature 
rake, 

wall temperatures, 
jet flow rates. 



jet nozzle 

✓ outer dia = 1 inch 


liquid 


Bulk 


7.03 in 


2.3 in 


5.75 in 


5 in 



NASA TM-106629, 1994 


pump section 
outer dia = 2.52 inch 


c\i 


K-site LH2 Axial Jet: Test Runs Simulated 




Table 1: Experimental conditions for mixing tests. 


Test run# 

Initial liquid fill, % 

Jet volume flow 
rate, m-Vhr 

Jet velocity, m/s 

mt 7 

Re Jet 

Initial pressure, kPa 

436 

85.3 

1.82 

1.32 

160,000 

187.0 

434 

86.3 

3.47 

2.51 

304.000 

186.1 

457 

49.1 

1.82 

1.33 

161,000 

186.5 

449 

49.1 

3.41 

2.47 

299,000 

186.1 


1.82 m 3 /hr = 8.0132 GPM LH2 
3.41 m 3 /hr = 15.0138 GPM LH2 
3.47 m 3 /hr = 15.278 GPM LH2 
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P(kPa) 
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K-site LH2 lg Axial Jet: Experiment and SINDA/FLUiNt 


SINDA/FLUINT used 75 LIQUID LUMPS, 50 VAPOR LUMPS 
Test Run 436 (85.3% liquid fill, 8.0132 GPM LH2 jet flow rate) 



M£SSU*£ VS TIME 
HEAT JUTE INTO TANK - 4.2 W,W*2 
K5TTE 1ET CASE 4* 



7.4 7.45 7.5 7.5S 7.5 7.65 7.7 7.75 7.B 7.B5 7.« 7.95 B B.05 8-X 8.15 8.2 8.25 8.3 8.3S B.4 

time (hr) 




P(kPa) 






K-site LH2 lg Axial Jet: Experiment and SINDA/FLUilMt 


'•InWsa 


SINDA/FLUINT used 75 LIQUID LUMPS, 50 VAPOR LUMPS 
Test Run 434 (86.3% liquid fill, 15.278 GPM LH2 jet flow rate) 


WES9J&E V5 TIME 
HEAT HATE I WTO TANK a 4.2 
KSfTE )£T CASE 4 >4 



7.165 7.19 7.195 7.2 7.205 7.2t 7.215 7.22 7-225 


7.2SS 7-26 7.26$ 7.27 7.275 7.28 72BS 7.29 7.295 


-• TEST data 
— 5F75LB0V 


time (hr) 


P(kPa) 




K-site LH2 lg Axial Jet: Experiment and SiNDA/FLUJNt 


SINDA/FLUINT used 75 LIQUID LUMPS, 50 VAPOR LUMPS 
Test Run 457 (49.1% liquid fill, 8.0132 GPM LH2 jet flow rate) 
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PRESSURE VS TJVt 
HEAT RATE INTO TANK H.2 
KSTTF JET CASE 457 


PRESSURE (KPa) 



• • TEST DATA 
— Sf .75LHOV 


P(kPa) 
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K-site LH2 lg Axial Jet: Experiment andSIIViDA/FUJINT 


SINDA/FLUINT used 75 LIQUID LUMPS, 50 VAPOR LUMPS 
Test Run 449 (49.1% liquid fill, 15.0138 GPM LH2 jet flow rate) 


♦ rx 

• NA&tf 

v. 1 ■ ■ 


PRESSURE V5 TIME 
HEAT HATE INTO TANK = 4.2 VV’M~ 2 
KSHE JET CASE 449 



6. IS 8-2 8 75 S3 £ 35 8.4 B.45 8.5 B.55 B.6 8.65 8.7 8.75 8.8 BBS SO 6 05 0 


time (hr) 


